%%
%作业3
%1
x = 1:2/50:3;

p = x.^5-9.*x.^4+31.*x.^3-51.*x.^2+40.*x-12;
plot(x,p,"Marker","o")
%%
%2
clc,clear
syms x y
ezplot((abs(x))^(3/4)+(abs(y))^(3/4)-1);
hold on
ezplot((x^2+y^2)^3-(x^2)*(y^3))
hold off
%%
%3
clc,clear
%syms a b
x = linspace(-3,3,33);
y = linspace(-3,3,33);
[a,b] = meshgrid(x,y);
z = 3.*(1-a).^2.*exp(-a.^2-(b+1).^2)-10.*(a/5-a.^3-b.^5).*exp(-a.^2-b.^2)-(1/3).*exp(-(a+1).^2-b.^2);
surf(a,b,z);
grid on
%%
%4
clc
clear

f = @(x,y,z) x.^2+y.^3+z.^4-1;
interval = [-2 2 -2 2 -2 2];
fimplicit3(f,interval)


